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We investigate the "twist" mode (rotation of the upper against the lower hemisphere) of a dilute 
atomic Fermi gas in a spherical trap. The normal and superfiuid phases are considered. The linear 
response to this external perturbation is calculated within the microscopic Hartree-Fock-Bogoliubov 
approach. In the normal phase the excitation spectrum is concentrated in a rather narrow peak 
very close to the trapping frequency. In the superfiuid phase the strength starts to be damped and 
fragmented and the collectivity of the mode is progressively lost when the temperature decreases. 
In the weak-pairing regime some reminiscence of the collective motion still exists, whereas in the 
strong-pairing regime the twist mode is completely washed out. The disappearance of the twist 
mode in the strong-pairing regime with decreasing temperature is interpreted in the framework of 
the two-fluid model. 

PACS numbers: 03.75.Ss,21.60. Jz 



I. INTRODUCTION 

The experimental and theoretical development of Bose- 
Einstein condensation of trapped bosonic atoms Q has 
also triggered the investigation of trapped atomic Fermi 
gases at very low temperature . One of the main goals 
in the research of these Fermi systems is to detect the 
transition from the normal to the superfiuid phase, as- 
sociated with the appearance of a macroscopic order pa- 
rameter of strongly correlated Cooper pairs below a cer- 
tain critical temperature Tc- In order to have an attrac- 
tive s-wave interaction which can provide the pairing cor- 
relations, the atoms must be trapped and cooled in two 
different hyperfine states as has been achieved in several 
recent experiments (see, e.g., Ref. Q). From theoretical 
side, the pairing problem of trapped fermions has been 
studied from different points of view 

Besides the ground-state properties, there is also in- 
terest in knowing the spectrum of collective excitations. 
As we stated above, ultracold atomic Fermi gases are as- 
sumed to become superfiuid below Tc, and it is therefore 
important to study low energy collective modes also in 
the superfiuid phase 0, ll| . Having different properties 
in the normal and superfiuid phases, these excitations 
can serve as experimental signals for superfluidity. For 
instance, the frequencies of breathing modes of trapped 
atomic Fermi gases measured in recent experiments 
give strong indications that the superfiuid phase has been 
reached. 

It is interesting to compare the situation of trapped 
fermionic atoms to that of atomic nuclei, which can also 
show a superfiuid behavior. Contrary to the nuclear case, 
the fact that the interaction in atomic gases is tunable 
experimentally allows to study the collective modes in 
different regimes. For dilute systems, the atom-atom in- 
teraction can be parametrized by a zero-range force pro- 
portional to the s-wave scattering length between atoms 



in two hyperfine states 0, |^] . By changing the applied 
magnetic field around a Feshbach resonance the s- 
wave scattering length can be modified. In one limiting 
regime, that of weak pairing, which is similar to the situ- 
ation in atomic nuclei, the pairing results only in a small 
perturbation to the response of the system to the external 
probe. In the other limit of strong pairing the response 
is dominated by the effects of superfluidity. 

Many of the collective excitations show features proper 
to Landau's zero sound modes in bulk Fermi liquids 
[m which for finite Fermi systems translate into modes 
analogous to those of an elastic body Since the 

trapped atomic Fermi gases contain a very large num- 
ber of atoms, the single-particle orbital angular momenta 
near the Fermi surface can also become very large. Con- 
sequently, important orbital effects such as excitations 
having angular momentum and parity — 1"*" and 
2^ will exist, which correspond to magnetic resonances 
of Ml or M2 type, respectively, in atomic nuclei. The 
2^ excitation is the so-called twist mode, in analogy to 
the quadrupole torsional vibration of an elastic sphere 
[mil. From a macroscopic point of view, the twist 
consists of a coherent counterrotation of the particles in 
the upper hemisphere against those in the lower hemi- 
sphere. For small amplitudes, it corresponds to a purely 
kinetic excitation without spatial distortion of the equi- 
librium shape. 

The twist mode has been studied in different Fermi 
systems. In nuclei, this mode has been analyzed from a 
semiclassical point of view within a fluid-dynamical de- 
scription 13] . From a quantum mechanical point of view, 
this mode has been studied so far only for magic nuclei 
(i.e., without pairing) such as ^''Zr and ^ospb [ll [l5l|. 
More recently, some experimental effort has been done to 
detect this mode by backward inelastic electron scatter- 
ing ^3 . A direct evidence for the existence of the orbital 
twist mode (to be distinguished from the 2~ spin-flip 
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mode) in nuclei has been achieved by comparing: electron 
and proton scattering cross-sections of ^*Ni The 
twist mode has also been theoretically studied in metal- 
lic clusters Ts^ although it has not yet been detected. 

So far, the theoretical study of the twist mode in 
trapped atomic Fermi gases has been done in the hy- 
drodynamical description and in the normal phase only 
P^] . In the case of a s-wave interaction, a moderate shift 
of the twist frequency of about 10% with respect to the 
non-interacting case was found, which is consistent with 
the fact that for a transverse zero-sound the s-wave inter- 
action does not contribute to the restoring force 0, . 

In the present article, our aim is different. We will an- 
alyze the effect of pairing correlations on the twist mode. 
This effect has not been considered in any of the theoret- 
ical studies mentioned above, neither for atomic nuclei 
nor for metallic clusters or trapped Fermi gases. Of par- 
ticular interest can be the study of the strong pairing 
regime, because it is known that in this case the low- 
energy collective modes are strongly affected by the pair- 
ing and can become signatures that the superfluid phase 
is reached 0, ll| ■ 

The paper is organized as follows. In Sec. 2, we 
sketch the derivation of the twist response function in 
the superfluid phase, using a Hartree-Fock-Bogoliubov 
or Bogoliubov-de Gennes framework [2lll2^. In Sec. 3, 
we consider the twist mode in the normal phase within 
a quantum-mechanical description. Sec. 4 is devoted to 
the study of the twist mode in the superfluid phase in 
the cases of weak and strong pairing correlations. Fi- 
nally, our conclusions are laid in the last section. 

II. QUASIPARTICLE RESPONSE FUNCTION 

In this article we will consider an atomic Fermi gas 
(atomic mass m), trapped in a spherical harmonic trap 
with frequency O. We assume that the atoms equally oc- 
cupy two hyperfine states, denoted by ct = ±1. Because 
of the low density of the gas, the interaction between 
the atoms can be regarded as pointlikc and its strength 
can be parametrized by the s-wave atom-atom scatter- 
ing length a. In order to simplify the notation, we will 
express all quantities in harmonic oscillator (h.o.) units, 
i.e., frequencies in units of fi, energies in units of Ml^ tem- 
peratures in units of Ml/ks, and lengths in units of the 
oscillator length Iho = \/fi/ {mVt). Furthermore, instead 
of the scattering length we will use the coupling constant 
g = Aira/lho as parameter of the interaction strength. 

The twist is a motion where the upper and lower 
hemispheres rotate in the opposite sense back and forth 



around the z axis with an angle proportional to z. This 
mode can be excited in both spherical and deformed 
(with a rotation axis) systems. Such a motion can be gen- 
erated by the operator zL^, where — —i{xVy — x) 
denotes the z component of the angular momentum op- 
erator. Restricting our description to small amplitudes, 
we can use linear response theory in order to treat the 
oscillations around equilibrium. Then the main problem 
consists in calculating the equilibrium state. In order 
to describe the system in the superfluid phase, this is 
done within the framework of a Hartree-Fock-Bogoliubov 
(HFB) or Bogoliubov-de Gennes [2lll23 | calculation sim- 
ilar to that presented in Ref. 0, but with the modified 
regularization scheme for the gap equation described in 
Refs. "23, "231. We refer to Ref. [23| for more details about 
our approach. The calculation provides the wave func- 
tions Ua{r) and Va{r) satisfying the HFB equations 

[Ho + W{r)]u^{r) + A(r)«„(r) = E^u^{r) , 
A(r)u„(r) - [i/o + W{r)]v^{r) = E^v^{r) . ^ ' 

Here Hq denotes the hamiltonian of the non-interacting 
h.o. minus the chemical potential, Hq ~ (— -|- 
r^)/2 — /i, while the interaction is accounted for in a 
self-consistent way through the Hartree potential W{r) 
and the pairing field A(r). 

Now let us consider the retarded correlation function 

poo 

Uo{u;)^-i cite^"*(([Q(t),Q(0)])), (2) 
Jo 

where ((•)) means the thermal average. In our case, Q is 
the twist operator 

QW=E I d?ri>l{t,r)zL,ij,{t,r). (3) 

The field operator -0 can be expressed in terms of quasi- 
particle creation and annihilation operators 6^ and h as 
follows: 

nlm 

-'^bl,n-.<Ur)e-^''-")- (4) 

Separating the radial and angular dependence of the 
wave functions, M„;m(r) = Uniir)YimiO, (/)) and Vnimir) = 
Vni{i")Yijn{0,(j)), one obtains after a straight-forward but 
tedious calculation the following result: 
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no(tj) = 2 51 ™^ / dnYi*^{9, (I)) cose Yi 
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iEn'i'~E,u)[f{Eni)~f{En'i')] 
{u; + irif-{E,,,i'~E,,iY 



drr'^[uni{r)un'i'{r) +Vni{r)vn'i'{r)]\ . (5) 



In deriving this formula, we have used the anticommuta- 
tion relations between the operators b and [{ba, bp} — 
— 0, {ba,b''ij} = Sa/i] as well as the relation 
{{bl^bp)) = f{Ea)Sai3, where / denotes the Fermi dis- 
tribution function, f{E) = l/(e^/^ + 1). Note that the 
relative signs appearing in the radial integrals in Eq. (O 
are different from those obtained, e.g., for the case when 
Q is a multipole operator as in Ref. @. The reason is 
that the twist operator is odd under time reversal, i.e., 
/ d''rrir)zL,g{r) = -[/ d^rg* ir)zL,f{r)]* . The angu- 
lar matrix element in Eq. jSJl can be computed explicitly, 
with the simple result 
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dnY*ie,c^) cos 0Yi,. 



otherwise . 



(6) 



Therefore the numerical task of calculating XIq reduces es- 
sentially to calculating the radial integrals and the triple 
sum over n, n' , and /. 

In general it is not sufhcient to calculate the free quasi- 
particle response IIq. Rather one has to calculate the 
QRPA response, which accounts for correlations with 
the quantum numbers corresponding to the excitation 
under consideration in the ground state. However, be- 
cause of the particular form of the interaction used here, 
it is clear that there cannot be any ground state cor- 
relations with the quantum numbers of the twist mode 
{J^ ~ 2^). Therefore the QRPA response function, 11, 
is just equal to the free quasiparticle response function. 
Ho 01 . In this sense the situation for trapped atoms is 
different from that in nuclei, where the spin-orbit part of 
the interaction leads to a (small) change of the twist re- 
sponse function e.g., through the coupling between 
the twist mode and the spin-flip mode, which is excited 
by the operator (r ® <t)2o- 

In the remaining part of this article we will show 
numerical results for the strength function S{u!) = 
— Im II{uj)/tt which we calculate from Eq. (O with a finite 
width rj for each peak. 



III. NORMAL PHASE 

Let us first look at the normal phase of the system at 
zero temperature, i.e., we artificially put A = in Eq. 
Q. For this situation, there exist microscopic descrip- 
tions of the twist mode in nuclei ITsj and in metal 
clusters However, for the twist mode in trapped 

atomic gases there exists only a calculation 0| following 
the fluid dynamical approach developped by Holzwarth 
for the nuclear case ^ISj . This fluid dynamical approach 
allows to predict the twist frequency, but it cannot an- 
swer the question if the twist mode as a collective motion 
exists at all .13.] . 

In the case of a non-interacting h.o., it is straight- 
forward to show that the operator zL^ excites only tran- 
sitions with Lu — 1 (in units of Ml). Therefore, in the 
non-interacting h.o., the total strength is concentrated 
at w = 1. If now the Hartree potential W is switched on, 
two effects are to be expected: 

a) The energy difference between neighboring shells 
becomes larger (smaller) in the case of an attractive (re- 
pulsive) interaction. Therefore, the twist frequency will 
be shifted upwards (downwards). This effect has been 
described quantitatively within the fluid dynamical ap- 
proach [l9l |. 

b) The degeneracy of states with different / is lifted, 
and we therefore expect a fragmentation of the strength 
of the twist mode into many particle-hole states corre- 
sponding to transitions n, Z — > n, l+l and n, I l—l 
[remember that for given quantum numbers n and Z, the 
number of h.o. quanta is 2{n — 1) + Z]. 

Both effects can be observed in Fig. ^ where we dis- 
play the strength function S{uj) of the twist mode as 
a function of the excitation energy for two systems of 
®Li with different numbers of atoms (scattering length 
a — — 2160ao [l^, where ao is the Bohr radius) in a trap 
with a frequency oi ~ 2tt x 144 Hz, corresponding to 
a coupling constant g — —0.4 in h.o. units. In order 
to show the fragmentation of the mode, we display the 
response function in a small energy interval (containing 
100% of the total strength) with a very high resolution 
[t] — 10^"*). Let us first look at the result correspond- 
ing to 10^ particles in the trap (5 x Iff* particles per 
spin state, chemical potential fi = 62.6). One can clearly 
see that the average frequency is higher than 1 and that 
the strength is fragmented into two series of peaks, corre- 
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FIG. 1: Strength function S{uj) (xlO^; in h.o. units) for the 
twist mode in a gas of lO'' (sohd line) and 2 x 10^ (dashed 
hne) trapped ^Li atoms at zero temperature without pairing 
(A = 0) as function of the frequency to (in units of hfl). 



spending to the two series of transitions mentioned above 
under b). With 2 x 10^ particles (10^ particles per spin 
state, /i — 78.0), the Hartree field is stronger and there- 
fore both effects, fragmentation and shift of the average 
frequency, are enhanced. 

The difference shown by the strength at low frequencies 
of the systems containing 10^ and 2 x 10^ particles is 
related to the different single-particle spectra of these two 
systems. In the case of 10^ particles, the Hartree field W 
breaks the accidental degeneracy of the non-interacting 
h.o. single-particle levels, but the different h.o. major 
shells are still separated. However, when the number of 
particles in the trap grows, the Hartree field becomes 
strong enough to mix different h.o. major shells. This 
leads more or less accidentally to the fact that in the case 
of 2 X 10^ particles the transition energies of the series 
n,l n -I- 1, / — 1 with 2{n — 1) + I = 82 (the major 
shell number 82 is the last one lying completely below 
the Fermi level) are almost degenerate at w 1.05. 

In order to compare our results quantitatively with 
the predictions obtained within the fluid-dynamical ap- 
proach, we define an average frequency according to 
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(7) 



For both cases considered here, this average frequency 
is in perfect agreement with the frequency Ufd predicted 
in Ref. 19] in the framework of the fluid-dynamical ap- 
proach: for N = 10^ atoms, uJav — 1.088 and Ufd = 1.087, 
and for TV = 2 X 10^ atoms, uJav = 1.100 and ujfd = 1.101. 
It should also be emphasized that the width of the inter- 
val over which the strength is distributed is very narrow 
compared with the average frequency of the twist mode. 



FIG. 2: Strength function S{u) {xlO''; in h.o. units) for the 
twist mode in a system with about 1800 atoms of ®Li at T = 
(top), T = 0.13 (middle) and T = 0.15 (bottom; lo and T in 
units of hfl) 



It is therefore justified to speak about a collective exci- 
tation. 



IV. SUPERFLUID PHASE 

Let us now consider the superfluid case. It is interest- 
ing to analyse how the properties of the twist mode are 
modified when pairing correlations are taken into account 
in the calculations and the full HFB equations are solved. 
We will show that the structure of the strength function 
and the collectivity of the twist mode are strongly af- 
fected by pairing correlations and we will study this de- 
pendence at different temperatures for two systems with 
different numbers of atoms. 

We set the coupling constant g equal to —1 in h.o. 
units. For ^Li atoms with a scattering length a = 
— 2160ao this corresponds to a trapping frequency of 
uj = 27T X 817 Hz. (We chose a stronger coupling than 
in the previous section in order to be able to study the 
case of strong pairing, which would be possible only for 
extremely large numbers of particles if <? = —0.4.) We 
shall consider two cases for the trapped gas: a) A small 
system with around 1800 atoms (weak pairing regime, 
A < hfl); b) A big system with around 3.6 x lO'' atoms 
(strong pairing regime, A > hTl). For both cases we will 
take into account different temperatures and analyze how 
the twist mode evolves when the critical temperature Tc 
of the phase transition is approached and crossed. 

Before passing to consider the two cases we would like 
to mention that in atomic nuclei, which are the only sys- 
tems for which the twist mode has been observed so far, 
one is always in the weak pairing regime, the relation 
A < hn being always satisfied. 

a) Weak pairing regime: The chemical potential /i in 
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this case is chosen equal to 16. We show in Fig. 13 the 
strength function for three values of the temperature: 
T = (top), T = 0.13 (middle) and T = 0.15 (bottom). 
The three cases correspond to values of the pairing field 
in the center of the trap of A(r = 0) = 0.67, 0.26 and 
0, respectively (in h.o. units). In the last case (bottom 
of the figure) the gas is in the normal phase: we ob- 
serve that in the normal phase the strength function is 
concentrated at about ui — 1.12 (this is slightly higher 
than in Fig. ^ because of the stronger coupling, which 
leads to a stronger Hartree field). If we lower the tem- 
perature, the supcrfluid transition takes place; the effect 
on the strength function is to push its structure towards 
higher values of the energy. Qualitatively this can be 
understood by replacing the single-particle energies e„; 
by the quasiparticle energies Eni « \/ {^ni — m)^ + 
where A denotes the average matrix element of the pair- 
ing field at the Fermi surface. Neglecting the effect of 
the Hartree field for the moment, one obtains in this 
way a shift of the twist frequency from 1 to a higher 
value which lies between Vl + and A + \/l + A^. 
To see this, let us consider two limiting cases: If the 
chemical potential lies exactly on a single-particle level 
(half- filled shell), /i = Np -1-3/2, a transition of the type 
Np — > Np + 1, for example, corresponds to the cre- 
ation of two quasiparticles with energies E^p = A and 
Emp+i = vT-j-A^. In the other limiting case, the chemi- 
cal potential lies between two single-particle levels (closed 
shell), /i = Np + 2, and the twist mode corresponds to 
the excitation of two quasiparticles having each the en- 
ergy Enp = En^+i = v/l/4 + A2. 

Moreover, as one can also observe in Fig. 13 the exci- 
tation mode becomes less collective and, due to pairing, 
more and more damped and fragmented if one goes from 
T — 0.13 to T = 0. In the latter case pairing correla- 
tions are more intense and the loss of collectivity and the 
Landau damping are consequently more important. A 
similar Landau damping effect due to superfluidity has 
been found in Ref. |3] for the spin-dipole mode in the 
weak pairing regime. 

An other interesting effect to notice is the strength 
below Lo = 1 which appears below Tc but disappears at 
T = 0. Obviously this effect is due to the second term in 
Eq. which is equal to zero at T = 0. 

b) Strong pairing regime: Let us consider now the case 
with about 3.6 x lO'' atoms (/i = 40). We present in Fig. 
O the strength function at four temperatures: T = 0, 
2, 5, and 6.5 (from top to bottom). In the four cases 
the central values of the pairing field are A(r = 0) = 
12.7, 12.5, 9.6, and 0, respectively. In the latter case 
(normal phase) we observe a unique peak centered at 
about Ld — 1.2. Again, the energy is higher with respect 
to Fig. n and with respect to the case a) (Fig. O due 
to the stronger coupling. A fragmented structure with a 
very low strength exists in the energy region from uj — 
2 up to = 4. When we lower the temperature, we 
cross the superfluid transition (see upper panel of Fig. 
0}. At r = 5 the system is superfluid: we observe that 
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FIG. 3: Strength function S{lu) (xlO^; in h.o. units) in a 
gas with about 3.6 x 10* atoms at four different temperatures 
(from top to bottom): T — 0, 2, 5, and 6.5 {uj and T in units 

of hn). 

the main peak still exists, even if the excitation mode 
is less collective than in the normal phase case. Also, 
the fragmented structure towards w = 2 — 3 gets more 
strength than in the previous case. The fragmentation 
becomes much stronger and extends up to w = 10 when 
the temperature is lowered further, as can be seen in the 
case T = 2. However, the peak at about uj — 1.2 is still 
visible. Finally, if we look at the T = case, where the 
pairing correlations are the strongest, we observe that 
the main collective peak completely disappears, while a 
very fragmented structure with a low strength remains 
in the energy region between uj = 2 and uj = 10. We can 
thus conclude that at T = the collective twist mode 
does not exist any more. The same conclusion has been 
drawn in Ref. 8J for the spin-dipole mode in the strong 
pairing regime. 

Actually, once the irrotational flow limit (strong pair- 
ing) is reached 26, 27], the superfluid current has an 
irrotational velocity held, and the only possible excita- 
tions of the superfluid are density-fluctuation modes. In 
the language of a two fluid model, all the other excita- 
tions of the gas, such as the twist and the spin-dipole 
modes, have to be related to its normal component, as 
was discussed in Refs. |2^|23. When one decreases the 
temperature below Tc, the number of "normal" quasipar- 
ticles is reduced and therefore the strength of the twist 
mode becomes smaller. On the other side, the energy 
spectrum of the normal quasiparticles is modifled, lead- 
ing to a destruction of coherence between quasiparticles 
moving in the same direction It follows that the 

mode is more and more damped when one approaches 
T = 0. Obviously, this effect will strongly depend on the 
strength of pairing correlations, and this is why it is more 
important in the strong pairing regime. 
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Let us now discuss the relationship between the 
strength of the twist response function and the normal- 
fluid component of the system in a more quantitative way. 
To that end we consider the inverse-energy weighted sum 
rule, which is proportional to the real part of the response 
f miction at cj = 0: 

r.„3i£).4n,o). (8, 

Jo 2 

Within the two-fluid model it can be shown explicitly (see 
appendix) that this quantity is related to the density of 
the normal-fluid component of the system, /9„, by 

In Fig. 0] (bottom) we show numerical results for the 
dependence of the sum rule on temperature. The solid 
line represents the full HFB calculation, while the dashed 
line corresponds to Eq. The agreement is very sat- 
isfactory except at extremely low temperature, where 
quantum finite-size effects (corrections in hu/A, see Ref. 
|27|) lead to a non- vanishing value of the sum rule, 
whereas the two-fluid model predicts that the sum rule 
should go to zero at zero temperature because of the 
vanishing normal-fluid component. However, the overall 
good agreement confirms our interpretation that only the 
normal-fluid part of the system participates in the twist 
motion. In order to recognize more easily the regions 
where the gas is superfluid and normal, and to observe 
how pairing correlations decrease by increasing the tem- 
perature, we also plot in Fig. 0] (top) the value of the gap 
at the center of the trap, A(0). Note that the temper- 
ature dependence of 11(0) differs considerably from that 
of A(0). 

To conclude this section, we remark that our approach 
is only valid in the regime where collisions between atoms 
can be neglected. Following Ref. , this means that the 
mean time between collisions, r = pa'^vp(T / tp)"^ (where 
p, vp, and ep are the density, Fermi velocity, and Fermi 
energy, respectively) , must be much larger than the oscil- 
lation period in the trapping potential, 27r/r2. Expressed 
in h.o. units, we obtain ilT/{2Tr) = 6TT^/{gT)'^. In the 
case of T = 6.5 this ratio gives still 4.4, i.e. an atom per- 
forms more than four oscillations before it collides with 
another one. Consequently, all cases we considered are 
well in the coUisionless regime. 

V. SUMMARY AND OUTLOOK 

In this article we have studied the twist mode of 
an atomic Fermi gas trapped by a spherical harmonic 
potential in the normal and in the superfluid phase. 
The ground state has been obtained by solving the 
Bogoliubov-de Gennes equations, usingthe regulariza- 
tion procedure introduced in Ref. [23l 12^ . The excita- 
tions have been treated within the linear response theory. 




FIG. 4: Temperature dependence of the gap in the center 
of the trap, A(0) (top; A and T in units of hfl), and of the 
static response function —11(0) (bottom; xlO~^inh.o. units), 
which is equal to twice the inverse-energy weighted sum rule, 
for a gas with 3.6 x 10* atoms. In the lower figure, we show 
for comparison the result of the HFB calculation (solid line) 
together with the result obtained within the two-fluid model, 
Eq. @ (dashed line). 

As the zero-range interaction does not couple to the twist 
operator, we analyzed this excitation by calculating the 
free quasiparticle response function. 

We have analyzed the twist mode without pairing cor- 
relations by setting A = in the Bogoliubov-de Gennes 
equations. We observed that the strength function is 
concentrated around an energy higher than uj = I. This 
shift (with respect to the case of a non-interacting h.o.) 
is due to the Hartree potential and depends on the sign 
of the coupling constant g. We have also observed a frag- 
mentation of the strength which describes the transitions 
n,l —^ n,l + 1 and n,l n + 1,1 — 1. 

In the case of pairing correlations we have shown that 
the excitation mode starts loosing its collectivity below 
the critical temperature Tc- When the temperature is 
lowered from Tc towards T = 0, the strength function 
becomes more and more damped and fragmented. In the 
weak pairing regime (A < hil) this effect is less pro- 
nounced than in the strong pairing regime (A > HQ): 
In the weak pairing case the collective twist mode still 
exists at zero temperature. With increasing strength of 
the pairing correlations, the collective peak is shifted to 
higher energies, and at the same time it becomes more 
and more broad and fragmented and its strength de- 
creases. Finally, in the strong pairing limit it completely 
disappears at T = 0. In fact, it can be predicted that 
the twist mode ceases to exist once the pairing is strong 
enough for the system to reach its irrotational flow limit 
PI 12^. 

It should be pointed out that, in the normal phase, the 
twist mode can only exist in the coUisionless regime, since 



7 



the restoring force for this collective oscillation comes 
entirely from the Fermi surface deformation |13|. This 
means that detecting the twist mode in the normal phase 
would be a signal that the system is in the coUisionless 
regime. This might be of importance since the evidence 
for the superfluidity obtained in recent experiments 9j 
relies on the assumption that the system is in the col- 
lisionless regime. The subsequent disappearance of the 
twist mode at lower temperatures would be a clear signal 
that the superfluid phase has been reached. Concerning 
the possibility to excite the twist mode experimentally 
we refer to Ref. [1^ . 

Recently, the twist mode has been measured in open- 
shell finite nuclei such as ^^Ni In the existing the- 
oretical studies of the twist mode in nuclei pairing cor- 
relations have not been taken into account, i.e., these 
studies are essentially restricted to closed-shell (magic) 
nuclei. Although nuclei are in the weak pairing regime, 
we think that a theoretical study of the twist mode in 
nuclei taking into account pairing correlations could be 
very interesting. Work in this direction is in progress. 
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APPENDIX A: RELATION BETWEEN THE SUM 
RULE AND THE NORMAL COMPONENT 

In this appendix we will briefly show how a relation- 
ship between 11(0) and the normal-fluid component of a 
system with strong pairing (A 3> 1 in h.o. units) can be 
established. A detailed discussion of some of the topics 
mentioned here can be found in Ref. l23 . 



As mentioned in Sec. 4, the inverse-energy weighted 
sum rule is proportional to the response 11(0) of the sys- 
tem to a static perturbation with a perturbation hamil- 
tonian Hi cx zL^. By taking the fi — > limit of the time- 
dependent HFB equations, one can derive equations sim- 
ilar to the Vlasov equation for the superfluid phase (in 
our case, of course, the time-dependence does not play 
any role). The resulting deviation of the Wigner func- 
tion p{r,p) from its equilibrium value reads 

^^(^'^) = (^)...(./^(^'^)' 

with 

E{r,p)^ V[/i(r,p)]2 + A2(r) , (A2) 

where h{r,p) and hi{r,p) denote the Wigner transforms 
of Ho + W{r) and Hi, respectively. Since 11(0) is deflned 
as the expectation value of zLz in the perturbed system, 
we can write 

"<°'-/1^(^)....,,<— 

(A3) 

Assuming spherical symmetry and a strongly peaked 
Fermi surface (i.e., A,T <C /i) it is straight-forward to 
derive Eq. 0, where 

.„W^.(r)/.«(-^«^)^^^^^ (A4) 

is the density of the normal-fluid component within the 
two-fluid model. 

Note that the temperature dependence of Eq. © is dif- 
ferent from that of the number of normal particles, since 
in Eq. © the factor weights very strongly the surface 
of the system, where A(r) is smaller and where conse- 
quently the normal- fluid fraction Pn/ p is higher than in 
the center of the trap. 
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